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I. INTRODUCTION 

Nonequilibrium transitions from stuck to flowing phases underlie the physics of a wide range 
of physical phenomena. In a first class of systems the onset of a stuck or frozen state occurs 
as a result of intrinsic dynamical constraints, due to interactions or crowding, and is usually 
referred to as jamming ^ | . Familiar examples are supercooled liquids that become a glasses upon 
lowering the temperature, colloidal suspensions that undergo a glass transition due to crowding 
upon increasing the density or the pressure, foams and granular materials that jam under shear, 
arrays of dislocations in solids that jam under an applied load. In a second class of systems the 
transition to a stuck state is due to external constraints, such as the coupling to quenched disorder 
(pinning centers from material defects in vortex lattices, optical traps in colloids, etc.), and is 
denoted as pinning [jj. Both classes of systems can be driven in and out of glassy states by tuning 
not only temperature, density or disorder strength, but also an applied external force. The external 
drive may be a shear stress in conventional glasses or simply a uniform applied force in systems 
with extrinsic quenched disorder, where even a uniform translation of the system relative to the 
fixed impurities represents a nontrivial perturbation. Vortex lattices in superconductors [jj] and 
charge density waves (CDWs) in metals [j| can be driven in and out of stuck glassy states by a 
uniform external current or electric field, respectively. As recognized recently in the context of 
jamming, the external drive plays a role not unlike that of temperature in driving the system to 
explore metastable configurations and should be included as an axis in a complete phase diagram. 

In this lectures I will focus on zero-temperature depinning transitions of interacting condensed 
matter systems that spontaneously order in periodic structures and are driven over quenched 
disorder. The prototype examples are vortex lattices in type-II superconductors P| and charge 
density waves in anisotropic metals j^j . Other examples include Wigner crystals of two dimensional 
electrons in a magnetic field moving under an applied voltage p, lattices of magnetic bubbles 



moving under an applied magnetic field gradient Q|, and many others. In general these systems 
form a lattice inside a solid matrix, provided by the superconducting or conducting material and 
are subject to pinning by random impurities. The statics of such disordered lattices have been 
studied extensively P] . One crucial feature that distinguishes the problem from that of disordered 
interfaces is that the pinning force experienced by the periodic structure is itself periodic, although 
with random amplitude and phase [3]. As a result, although disorder always destroys true long- 
range translational order and yields glassy phases with many metastable states and diverging energy 
barriers between these states, the precise nature of the glassy state depends crucially on disorder 
strength. At weak disorder the system, although glassy, retains topological order (the resulting 
phase has been named Bragg glass in the context of vortex lattices) |9(. Topological defects 
proliferate only above some characteristic disorder strength, where a topologically disordered glass 
is formed. 

The driven dynamics of disordered periodic structures have been studied extensively by modeling 
the system as an overdamped elastic medium that can be deformed by disorder, but is not allowed 
to tear, that is by neglecting the possible formation of topological defects due to the competition 
of elasticity, disorder and drive. This model, first studied in the context of charge density waves, 
exhibits a nonequilibrium phase transition from a pinned to a unique sliding state at a critical 
value Ft of the driving force. This nonequilibrium transition displays universal critical behavior 
as in equilibrium continuous transitions, with the medium's mean velocity v acting as an order 
parameter P, 0, While the overdamped elastic medium model may seem adequate to describe 
the dynamics of driven Bragg glasses, many experiments and simulations of driven systems have 
shown clearly that topological defects proliferate in the driven medium even for moderate disorder 
strengths [ill . I12I lid . I14I lo^ . The dynamics near depinning becomes spatially and temporally 
inhomogeneous, with coexisting moving and pinned degrees of freedom. This regime has been 
referred to as plastic flow and may be associated with memory effects and even hysteresis in the 
macroscopic response. 

The goal of the present lectures is to describe coarse-grained models of driven extended systems 
that can lead to history-dependent dynamics. Such models can be grouped in two classes. In the 
first class the displacement of the driven medium from some undeformed reference configuration 



remains single- valued, as appropriate for systems without topo 
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ogical_ defects, but the interactions 
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IjJ, |20] ■ In the second class 



are modified to incorporate non-elastic restoring forces 
of models top olog ical defects are explicitly allowed by removing the constraint of single-valued 
displacements js ll . I22I l^ . Here we will focus on the first class and specifically consider driven 
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periodic media with a linear stress-strain relation, where the stress transfer between displacements 
of different parts of the medium is nonmonotonic in time and describes viscous-type slip of neigh- 
boring degrees of freedom. A general model of this t ype that encompasses many of the models 
discussed in the literature was proposed recently by us 18, ^(JQ]- Here slips between neighboring 



degrees of freedom are described as viscous force, that allows a moving portion of the medium to 
overshoot a static configuration before relaxing back to it. It is shown below that such viscous 
coupling can be considered an effective way of incorporating the presence of topological defects 
in the driven medium. Related models have also been used to incorporate the effect of inertia or 
elastic stress overshoot in crack propagation in solids jlfll . . The precise connection between the 
two classes of models has been discussed in Ref. j^ . 

In Section [U] we review the simplest example of depinning transition, obtained when non- 
interacting particles are driven through a periodic pinning potential. By contrasting the case of 
periodic and non-periodic pinning, we stress that care must be used in the definition of the mean 
velocity of the system. In Section IIII1 we first describe the generic coarse-grained model of a 
driven elastic medium that exhibits a continuous depinning transition as a function of the driving 
force from a static to a unique sliding state. Next we introduce an anisotropic visco-elastic model 
as a generic model of a periodic system driven through strong disorder. The model considers 
coarse-grained degrees of freedom that can slip relative to each other in the directions transverse 
to the mean motion, due to the presence of small scale defects (phase slips, dislocations, grain 
boundaries) at their boundaries, but remain elastically coupled in the longitudinal directions. The 
slip interactions are modeled as viscous couplings and a detailed physical motivation for this choice 
is given m section IfflCl Most of our current results for these type of models are for the mean-field 
limit and are presented in Section ITVl The studies carried out so far for finite-range interactions 
suggest that the mean-field theory described here may give the correct topology for the phase 
diagram, although there will of course be corrections to the critical behavior in finite dimensions 
[2^ . Finally, we conclude in Section by discussing the relation to other models described in the 
literature and the connection to experiments. 



II. DEPINNING OF NONINTERACTING PARTICLES 



It is instructive to begin with the problem of a single particle driven through a periodic pinning 
potential as the simplest illustration of driven depinning. Assuming overdamped dynamics, the 
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equation of motion for the position x of the particle is 

(^ = F + hY(x), (1) 

where £ is a friction coefficient (in the following we choose our units of time so that ( = 1), F 
is the external drive and Y(x) = Y(x + n), with n an integer, is a periodic function of period 1. 
For simplicity we choose a piecewise linear pinning force, corresponding to Y(x) = (1/2 — x), for 
< x < 1. In this case a periodic solution of Eq. (^Q) is obtained immediately in terms of the 
time T needed to traverse a potential well, or period. Introducing an arbitrary time tj such that 
if x(tj) = n, then x(tj + T) = n + 1, the particle position for tj + nT < t < tj + (n + 1)T is 

1 _ e -h(t-tj-nT) 

*(*)=" + f3 ^T . ( 2 ) 



where T is given by 



for i 7 > /i/2 and diverges for F < fa/2. In other words if F < h/2 the particle never leaves the 
initial well, i.e., it is pinned. The threshold force for depinning is then F c = h/2. In the sliding 
state the mean velocity is defined as the average of the instantaneous velocity v (t) = ^ over the 
arbitrary initial time t j. This gives 

v - (v) tj = d 4f v(t) = 1 . (4) 

Jt-(n+l)T 1 1 

This definition naturally identifies the mean velocity of the particle with the inverse of the period. 
The logarithmic behavior of v near threshold, v ~ —l/ln(F c — F), is peculiar to a discontinuous 
pinning force. For an arbitrary pinning force Y(X) the period T is 

f 1 1 

T = / dx — , (5) 

J F + hY(x) K > 

and can be evaluated analytically for various forces. For instance, for a sinusoidal pinning force, 
Y{x) = sin(27rx), one finds T = (F 2 — /i 2 ) -1 / 2 , which gives v ~ (F — F c ) 1 / 2 near threshold, a 
generic behavior for continuous pinning forces. 

The main focus of the remainder of this paper will be on the modeling of extended driven 
systems as collections of interacting degrees of freedom. It will then be important to distinguish 
two cases. For extended systems that are periodic, such as charge density waves and vortex 
lattices, the pinning potential is itself periodic as each degrees of freedom sees the same disorder 
after advancing one lattice constant. For non-periodic systems, such as interfaces, each degree of 
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freedom moves through a random array of defects. When interactions are neglected, an extended 
periodic system moving through a periodic random pinning potential can be modeled as a collection 
of iV non-interacting particles, where each particle sees its own periodic pinning potential. The 
pinning potentials seen by different particles may differ in height and be randomly shifted relative 
to each other, as sketched in Fig. 1. The equation of motion for the i-th particle at position X{ is 
then 

^- = F + h l Y(x l + 7i ), (6) 

where ji are random phases uniformly distributed in [0, 1) and the pinning strengths hi are drawn 
independently from a distribution p{hi). Since the displacements Xi are decoupled, they can be 
indexed by their disorder parameters 7 and h instead of their spatial label i, i.e., 7^, /ij) — ► 
x(t;j,h). The mean velocity of the many-particle system can then be written as an average over 
the random phases and pinning strengths, 

% 

= J dhp(h) J djv(t; 7 ,h) , (7) 

where v(t;^f,h) = • The average over the random phase of each degree of freedom is 

equivalent to the average over the random time shift tj described for the single-particle case and 
yields dj v(t; 7, h) = 1/T(h), with T(h) the period of each particle given in Eq. Q- The mean 
velocity is then 



> ( 8 ) 

T ( h )/h 

where {...)h = J dh...p(h) denotes the average over the barrier height distribution. For distributions 
p(h) that have support at h = 0, a system of noninteracting particles with periodic pinning depins 
at F = 0, as there are always some particles experiencing zero pinning force. 

A different single-particle problem that has been discussed in the literature is that of a particle 
moving through a random (non-periodic) array of defects 28(. The defects can be described as 
pinning potential wells centered at random positions and/or with random well heights. To make 
contact with the periodic case we consider a particle moving through a succession of evenly spaced 
pinning potential wells of random heights. The equation of motion is 

N p -1 

d t x = F+Y J h p Y{x-p), (9) 

p=0 



FIG. 1: (a) Sketch of nonintcracting degrees of freedom driven over a random periodic pinning potential 
in one dimension. Spatial coordinates have been discretized so that degrees of freedom are labelled by an 
index i. In (b) the case where each degree of freedom interacts elastically with its neighbors is shown. This 
is a discretized one-dimensional realization of the elastic medium model described by Eq. (|13J) below. 

where N p is the total number of pinning centers and the pinning strengths h p are drawn indepen- 
dently from a distribution p(h p ). Choosing again the piecewise-linear pinning force, the time to 
traverse the p-th well is simply T(h p ), with T given by Eq. (j3J). The mean velocity of the particle 
is defined as the total distance travelled divided by the total time and is given by 

y= ^ =^— . (10) 
E p T(h P ) (T(h)) h 1 j 

In this case, unless the distribution p(h) is bounded from above, there is always a finite probability 
that the particle will encounter a sufficiently deep potential well to get pinned. Therefore for 
unbounded p{h) the particle is always pinned in the thermodynamic limit. If p(h) is bounded from 
above by a maximum pinning strength h max , this value also represents the depinning threshold. 
Finally, the case of many noninteracting particles driven through a random array of defects is 
equivalent to that of a single particle, as the mean velocity of each particle can be calculated 
independently. The mean velocity of the system is then again given by Eq. (|lUj) . 




III. DEPINNING OF AN EXTENDED MEDIUM 



We consider a (i-dimensional periodic structure driven along one of its symmetry directions, 
chosen as the x direction. The continuum equations for such a driven lattice within the elastic 
approximation were derived by various authors by a rigorous coarse-graining procedure of the 



microscopic dynamics 
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311 ] . Assuming overdamped microscopic dynamics, the equation for 



the local deformation u(r, t) of the medium (in the laboratory frame) from an undeformed reference 
state is written by balancing all the forces acting on each portion of the system as 



d t Ui = djOij + F5 ix + F pi (r, u) , 



(11) 
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where is the stress tensor due to interactions among neighboring degrees of freedom, F is the 
driving force and F p is the periodic pinning force. The periodicity of the pinning force, which 
contains Fourier components at all the reciprocal lattice vectors of the lattice, arises from the 
coupling to the density of the driven lattice. 



A. Elastic model 



For conventional short-ranged elasticity the stress tensor is 

afj = 2cmUij + <%(cn - cm)ukk , (12) 

where en and cqq are the compressional and shear moduli of the driven lattice, respectively, and 
Uij = 7}(diUj + djUi) is the strain tensor. It was shown in Ref. [3] that deformations of the driven 
lattice along the direction of the driving force grow without bound due to large transverse shear 
stresses that generate unbounded strains responsible for dislocation unbinding. For this reason, 
we focus here on the dynamics of a scalar field u x (x,y,t) = u(r,t), with r = (x,y), describing 
deformations of the driven lattice along the direction of mean motion. The d— 1-dimensional vector 
y denotes the coordinates transverse to the direction of motion. Assuming en >> cqq, we obtain 
a scalar model for the driven elastic medium, given by 

d t u = c\\d 2 x u + c 66 VyU + F + F p (r, u) , (13) 

where F p denotes the x component of the pinning force. For simplicity we also consider a model 
that only retains the component of the pinning force at the smallest reciprocal lattice vector and 
choose our units of lengths so that the corresponding period is 1. The pinning force is then taken 
of the form 

F p (t,u) =/i(r)y(u(r,t)-7(r)) , (14) 

where Y(u) = Y(u + n) is a periodic function. The random pinning strengths h are drawn inde- 
pendently at every spatial point from a distribution with zero mean and short-ranged correlations 
to be prescribed below. The random phases 7 are spatially uncorrelated and distributed uniformly 
in [0,1). 

The model of a driven overdamped elastic medium embodied by Eq. (|13jl has been studied 
extensively both analytically and numerically P,Q, 33, It exhibits a depinning transition 



at a critical value Ft of the applied force from a static to a unique sliding state |35| . The depinning 
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can be described as a continuous equilibrium transition, with the mean velocity v = {dtu) playing 
the role of the order parameter, and universal critical behavior. The velocity vanishes as Ft is 
approached from above as v ~ (F — Fpy. The critical exponent (5 depends only on the system 
dimensionality and was found to be (3 = 1 — e/6+0(e 2 ) using a functional RG expansion in e = 4— d 

0,0. 

B. Viscoelastic model 

Strong disorder can yield topological defects in the driven lattice, making the elastic model 
inapplicable 0,0. In this case the dynamics becomes inhomogeneous, with coexisting pinned 
39 , |3| ■ The depinning transition may be discontinuous (first order) , possi- 



and moving regions 



bly with macroscopic hysteresis. Several mean-field models of driven extended systems have been 



proposed 



P,H,0,Q,lis 
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22 ] to describe this inhomogeneous dynamics. Here we focus on a 
class of models that retains a single- valued displacement field and a linear stress-strain relation, 
but assumes that the presence of topological defects can be effectively incorporated at large scales 
by a non-instantaneous stress transfer that couples to gradients of the local velocity (rather than 
displacement). More precisely, we consider an anisotropic model of coarse-grained degrees of free- 
dom that can slip relative to each other in at least one of the directions transverse to the mean 
motion, due to the presence of small scale defects (phase slips, dislocations, grai n boundaries) at 
their boundaries, but remain elastically coupled in the longitudinal directions [20|. This model in- 
corporates the anisotropy of the sliding state in the plastic flow region that results either from flow 
along coupled channels oriented in the direction of the drive (e.g., as in the moving smectic phase 
or in layered materials such as the high-T c cuprate superconductors. It also encompasses 
several of the models discussed in the literature. 

For generality, consider a d = dn + (i^-dimensional medium composed of degrees of freedom 
that are coupled elastically in d| direction and can slip relative to each other in the remaining dj_ 
directions. The axis x along which the driving force is applied is along one of the d\\ directions. 
The equation of motion for the displacement u(r^,r±,t) is given by 

d t u = KV\u + r]V\v + F + F p (t,u);, (15) 

with v = dtu the local velocity. This model will be referred to as the visco-elastic (VE) model as 
it incorporates elastic couplings of strength K in d| directions and viscous couplings of strength 
controlled by a shear viscosity 77 in the remaining d± directions. A two-dimensional cartoon of this 
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FIG. 2: A two-dimensional realization of the anisotropic driven medium described in the text. Spatial 
coordinates have been discrctized in the figure so that degrees of freedom are labelled by indices t and i, 
respectively transverse and longitudinal to the direction of the driving force, F. Each degree of freedom 
interacts with its neighbors via elastic couplings in the longitudinal direction and via viscous or similar slip 
couplings in the transverse direction. 

anisotropic model is shown in Fig. |2J 

For 77 = (or d± = 0) the VE model reduces to the elastic model (but with isotropic elasticity) 
of Eq. Conversely, for K = (or d\\ = 0) Eq. (|15j) reduces to the purely viscous model 

studied earlier by us [la . I2J]. For any distribution of pinning strengths with support at h = 0, 
the purely viscous model has zero threshold for depinning, but it does exhibit a critical point 
separating regions of unique and multivalued solutions for the mean velocity. In the VE model 
(77 7^ and K 7^ 0) even when fluid-like shear takes place, particle conservation gives a sharp 
depinning transition in flow along the channels. Furthermore, as shown below, the model has a 
sharp mean-field tricritical point separating a region of parameters where depinning is continuous, 
in the universality class of elastic depinning, from one where depinning become discontinuous and 
hysteretic. 

It is important to stress that the VE model still assumes overdamped microscopic dynamics. 
Velocity or viscous couplings can appear generically in the large-scale equations of motion upon 
coarse-graining the microscopic dynamics of a dissipative medium. In fact, next we show that 
viscous couplings indeed represent an effective way of incorporating the local dissipation due to 
the presence of topological defects. 



C. Viscoelastic coupling as an effective description of topological defects 

The goal of this section is to provide some justification to the anisotropic VE model as an 
effective description of topological defects in a driven lattice. To this purpose we consider a two 
dimensional medium and take advantage of the continuum equations developed many years ago 
by Zippelius et al. [^I to describe the time-dependent properties of two-dimensional solids near 
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melting. These authors combined the equations of free dislocation motion with solid hydrodynamics 
to construct a semimicroscopic dynamical model of a solid with free dislocations. They further 
showed that the dynamics of such of a "heavily dislocated solid" (an elastic medium with an 
equilibrium concentration of free dislocations) is identical to that of the hexatic phase obtained 
when a two-dimensional solid melts via the unbinding of dislocations More recently we [4J] 
reconsidered the dynamical equations for the "heavily dislocated solid" of Ref. [3] and showed 
that they can be recast in the form of the phenomenological equations of a viscoelastic fluid (with 
hexatic order) introduced many years ago by Maxwell j^ . In the presence of free dislocations 
the local stresses in the medium have contributions from both elastic stresses and defect motion. 
The latter couple again to the the local strains which control the defect dynamics. By eliminating 
the defect degrees of freedom, one obtains a linear, although nonlocal, relation between strain and 
stress, given by Q 

<r^ E (r, t) = Sij c L u kk (r, t) + Sijicn - c L ) [ dt'e"^^ v kk {r, t') 

+2c 66 f cft / e-(*-*')/ T '[ui J -(r,0-^«Ufcfc(r,0], (16) 
J -00 ^ 

where V{j = ^{diVj + djVi and the velocity v is defined here in terms of the momentum density 
g as v = g//9oi with po the equilibrium mass density of the medium. Also in Eq. ()16|) cl is 
the compressional modulus of the liquid and t\, « (cn/i^n/Oo)" 1 and t s rj (c66^ n / a o) _1 are 
the compressional and shear relaxation times, with /^' c the dislocation glide and climb mobility, 
respectively. Of course in the presence of dislocations the displacement u is no longer single- 
valued (although the strain Uij remains single-valued and continuous) and dtu 7^ v due to both the 
motion of vacancy/interstitial defects and of dislocations. The phenomenological Maxwell model 
of viscoelasticity is obtained by assuming that dtu = v even in the presence of dislocations. Then 
for t << t s , Tfe the viscoelastic stress cr VE (r,t) reduces to the familiar elastic stress tensor given in 
Eq. (H2), 

a YE (r,t«r s ,n)^a!}. (17) 

Conversely for t > > r s , r& one obtains 

<T^ E (r, t » T s ,n) w Sij CLU kk + 6ij(rjb + r])v kk + 2-qvij , (18) 

which describes stresses in a viscous fluid of shear viscosity 77 = cqqt s and bulk viscosity rjb = 
(en — cl)tij. The first term on the right hand side of Eq. (|18|) is the pressure and incorporates the 
fact that even a liquid has a nonzero long-wavelength compressional elasticity, which is associated 
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with density conservation. As we will see below this terms plays a crucial role in controlling the 
physics of depinning of a viscoelastic medium. The Maxwell viscoelastic fluid has solid-like shear 
rigidity at high frequency, but flows like a fluid at low frequency. Since the relaxation times t s and 
Tfe are inversely proportional to the density n$ of free dislocations, the Maxwell model behaves as 
a continuum elastic medium on all time scales when nj — > and as a viscous fluid when n^a\ ~ 1. 

Dislocation climb is much slower than dislocation glide (/iS << A^), resulting in r& >> r s . 
We therefore assume that the response to compressional deformations is instantaneous on all time 
scales, but retain a viscoelastic response to shear deformations. Letting — > oo, we find 

^ E (r, t) « Sij c lx u kk {t) + 2c 66 f d£e-W T * [ VlJ (t') - ^-v kk (t')] . (19) 

We now turn to the case of interest here, where topological defects are generated in a an 
extended medium driven through quenched disorder. In this case the medium has no low frequency 
shear modulus, but particle conservation still requires long wavelength elastic restoring forces to 
compressional deformations. On the other hand, the number of topological defects is not fixed 
as dislocations are continuously generated and annihilated by the interplay of elasticity, disorder 
and drive 0, Q, Furthermore, unbound dislocations can be pinned by disorder and do 
not equilibrate with the lattice. In the plastic region near depinning the dynamics remains very 
inhomogeneous and fluid-like and the pinning of dislocations by quenched disorder is not sufficient 
to restore the long wavelength shear-stiffness of the medium. For this reason we propose to describe 
the effect of topological defects near depinning by replacing elastic shear stresses by viscoelastic 
ones, while retaining elastic compressional forces. Of course the resulting model that assumes a 
fixed density of dislocations becomes inapplicable at large driving forces where dislocations heal as 
the lattice reorders. For the case of interest here of a scalar model describing only deformations 
along the direction of motion, the viscoelastic model of a driven disordered medium is 

d t u = c lx d 2 x u + c 66 [ dt'e- {t - t,)/T °d%v{t') + F + h(r)Y(u - 7 (r)) , (20) 

J — oo 

with v = dtu. This model naturally incorporates the anisotropy and channel-like structure of 
the driven medium, where shear deformations due to gradients in the displacement in the direc- 
tions transverse to the mean motion (d y u ^ 0) are most effective at generating the large stresses 
responsible for the unbinding of topological defects. It is instructive to note that due to the ex- 
ponential form of stress relaxation the integro-differential equation ()2()D is equivalent to a second 
order differential equation for the displacement, 

r s d\u + leS d t u = cudlu + ijajw + F + h{r)Y(u - 7 (r)) , (21) 
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with 7 c g an effective friction 



lgj ]. In other words the effect of a finite density of dislocations in 



the driven lattice yields "inertial effects" on a scale controlled by the time t s ~ 1/fi/- The purely 
viscous model obtained from Eq. (|21|) with en = was analyzed in detail in Ref. 2^1 where 
it was shown that if r s and r/ = cqqt s are tuned independently, then t s is a strongly irrelevant 
parameter in the RG sense. This allows us to consider a simplified form of the equation for the 
driven medium obtained from Eq. (|21|) with r s = 0, but r/ = cqqt s finite, leading to the general 
anisotropic viscoelastic model introduced in Eq. (|15|). 



IV. MEAN-FIELD SOLUTION 



The mean-field approximation for the VE model is obtained in the limit of infinite-range elastic 
and viscous interactions. To set up the mean field theory, it is convenient to discretize space in 
both the transverse and longitudinal directions, using integer vectors i for the du -dimensional intra- 
layer index and I for the (^-dimensional layer index. The local displacement along the direction 
of motion is u\{t). Its dynamics is governed by the equation, 

d t 4 = J2 - 4) + E w«i - + F + w(4 - it) , (22) 

0') < m > 

where the dot denotes a time derivative and (j) ((m)) ranges over sites j (m) that are nearest 
neighbor to % (I). The random pinning strengths h\ are chosen independently with probability 
distribution p(h % ^) and the 7^ are distributed uniformly and independently in [0, 1). For a system 
with N = iVy x N± sites, one mean field theory is obtained by assuming that all sites are cou- 
pled with uniform strength, both within a given channel and with other channels. Each discrete 
displacement then couples to all others only through the mean velocity, v = N^ 1 ^u\, and the 
mean displacement, u = N" 1 ^ u\. We assume that the disorder is isotropic and the system is self 
averaging and look for solutions moving with stationary velocity: u = vt. Since all displacements 
u are coupled, they can now be indexed by their disorder parameters 7 and h, rather than the 
spatial indices i and I. The mean-field dynamics is governed by the equation 

(1 + 77)12(7, h) = K(vt - u) +F + r]v + hY(u-j). (23) 

The cases K = and need to be discussed separately. 
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FIG. 3: (a) Velocity versus driving force for the purely viscous model (K = 0, rj ^ 0) with a narrow 
distribution of pinning strength, p(h) = 8{h — 1), for ?/ = 0,2,5. There is a finite dcpinning threshold at 
Ft = 1/2. In (b) the dcpinning and rcpinning forces and F± are shown as functions of rj. 



When K = 0, the mean field equation becomes identical to that of a single particle discussed 
in Section [H] driven by an effective force F + rjv (with friction 1 + rj) . In this case different degrees 
of freedom move at different velocities even in the mean field limit. The mean field velocity is 
determined by the self-consistency condition v = (u)^^, where the average over the random phases 
is equivalent to the average over the random times shifts tj given in Eq. @. For the case of a 
piecewise linear pinning force using Eq. (JSJ) we find 



with T(h,F) given by Eq. (J3J). The mean velocity obtained by self-consistent solution of Eq. (j24() 
is shown in Figs. and 0] for two distributions of pinning strengths. 

For a narrow distribution, p(h) = 5(h — 1), there is a finite threshold Ft = 1/2, independent of 
?7. The velocity is multivalued for any finite rj. When the force is ramped up adiabatically from 
the static state the system depins at Ft = Ft- When the force is ramped down from the sliding 
state, the system repins at the lower value F^(rj). The depinning and repinning forces are shown 
in Fig. O^b). The region where unique and multivalued velocity solutions coexist extend to rj = 0. 
For a broad distribution with support at h = 0, e.g., p(h) = e~ h , the threshold for depinning is 
zero as some of the degrees of freedom always experience zero pinning and start moving as soon as 
a force is applied. There is a critical point at (F c , rj c ). For rj > rj c the analytical solution for v(F) is 
multivalued, as shown in Fig. 0J If the force is ramped up adiabatically from zero at a fixed i] > rj c , 
the system depins discontinuously at i*| (77) , while when the force is ramped down it repins at the 
lower value (rj) , as shown in Fig. ^ The viscous model has also been studied in finite dimensions 



A. Mean-field theory for viscous model: K = 0, rj 7^ 




(24) 
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FIG. 4: (a) Velocity versus driving force for the purely viscous model (K = 0, r\ ^ 0) with a broad 
distribution of pinning strength, p(h) = e~ h for rj — 6, 16. In this case there are no stable static (pinned 
states). The velocity is single valued for r\ <r\ c and multi- valued for rj > r\ c . In this case when F is ramped 
up from zero, the velocity jumps discontinuously at where the system goes from the " slow- moving" to the 
"fast-moving" state. Here and below "coexistence" refers to multistability of the solutions to the equations 
of motion. When F is then ramped down from within the fast-moving state the jump in v occurs at the 
lower value F±. The forces F± and F^ become equal at the critical point, as shown in frame (b). 



by mapping it onto the nonequilibrium random field Ising model (RFIM) 2J| . In the mapping, the 
local velocities correspond to spin degrees of freedom, the driving force is the applied magnetic field 
and the mean velocity maps onto the magnetization. The RFIM has a critical point separating a 
region where the magnetization versus applied field curve displays hysteresis with a discontinuous 
jump to a region where there is no jump in the hysteresis curve 4|| 13- I n the viscous model 
the critical point separates a region where the velocity curve is smooth and continuous from the 
region where the "depinning" (from "slow-moving" to "fast-moving" states) is discontinuous and 
hysteretic. The critical point is in the Ising universality class, with an upper critical dimension 
d c = 6. 



B. Mean-field theory for VE model: K ^ and rj ^ 

When K ^ 0, all degrees of freedom are coupled by a spring-like interaction (the first term on 
the right hand side of Eq. I|23|)) to the mean field u = vt and cannot lag much behind each other. 
This forces all the periods to be the same, independent of h, and yields a nonvanishing threshold 
for depinning. In this case the mean field velocity is determined by imposing {u(t; 7, h) — vt)™^ = 0. 

It is useful to first review the case where i(/0 and rj = 0. In this limit, Eq. (|23|) reduces to 
the mean field theory of a driven elastic medium worked out by Fisher and collaborators ^(|. No 
moving solution exists above a finite threshold force Ft- For the piecewise linear pinning force this 
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is given by 

* = <5(FTSj>»- (25) 

For F > Ft there is a unique moving solution that has a universal dependence on F near i<y, 
where it vanishes as v ~ (F — Ft)^ ■ In mean- field the critical exponent j3 depends on the shape 
of the pinning force: (3 = 1 for the discontinuous piecewise linear force and [3 = 3/2 for generic 



1(| showed that 



smooth forces. Using a functional RG expansion in e = 4 — d, Narayan and Fisher 
the discontinuous force captures a crucial intrinsic discontinuity of the large scale, low-frequency 
dynamics, giving the gener al result P = 1 - e/6 + C(e 2 ), in reasonable agreement with numerical 



simulations in two and three dimensions 



33j, |34j. For simplicity and to reflect the "jerkiness" of 



the motion in finite-dimensional systems at low velocities, we use piecewise linear pinning below. 

When n > the nature of the depinning differs qualitatively from the rj = case, in that 
hysteresis in the dynamics can take place. Again, no self-consistent moving solution exists for 
F < Ft, with Ft independent of n. Above threshold, both unique and multi- valued moving 
solutions exist, depending on the values of the parameters: n, K, and the shape of the disorder 
distribution, p(h). To obtain the mean field solution in the sliding state, we examine the motion 
during one period T = 1/v during which the displacement advances by 1. Eq. (|23|) is is easily 
solved for < u < 1 and 7 = 0, with the result, 

Kvt + F + rjv + h/2 Kv . Xf 

= - h) = (i + ")A -JTTW + Ae ' (26) 

where A = {K + h)/(l + rj). At long times, regardless of the initial condition, u(t) approaches a 
periodic function of period T = 1/v with jumps in its time derivative at times tj + nT, with n an 
integer. The constant A is determined by requiring that if u(t j + nT) = n, then u(t j + (n + 1)T) = 
n + 1. Writing u(t; 7, h) = vt + u, it is easy to see that for an arbitrary value of 7, the solution u 
will have the form u = uivt — 7, h). The mean velocity is then obtained from (u(vt — 7, h))jh = 0. 
Averaging u over 7 is equivalent to averaging u for a fixed 7 over a time period, T, with the result, 



rtj + (n+l)T df 

(«>7 = / 777 u(v-jt,h) 

Jtj+nT 1 

F + rjv h 2 (K + 2h)v h 2 1 



K 2K(K + h) A (if + /i) X(if + /i) e A / I? - 1 
Finally, averaging over /i and using the consistency condition, we obtain 



(27) 



F — Ft = v[l — M(„, if)] + ( ^' + e(g+t)/( ; + „, _ 1 \ , (28) 
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FIG. 5: (a) Velocity versus driving force for the VE model with K = 1 and a broad distribution of pinning 



strength, p(h) 



\ The velocity is continuous and single- valued for r\ < n c and becomes multivalued for 



n > ?y c . The dashed line on he curve for 77 = 15 indicates the value F± where the system repins when the 
drive is ramped down from the sliding state. Frame (b) shows the depinning and repinning forces F^ and F± 
as functions of n. The tricritical point at (F c ,n c ) separates continuous from hysteretic depinning. Pinned 
and sliding states coexist in the region F± < F < F^. 



with Ft the threshold force given in Eq. (|25|) and M given by 

h 2 \ 

W+W/h 



M(rj,K) = (l + v) 



(29) 



As in the purely elastic case (77 = 0) only static solutions exist for F < Ft- For F > Ft there 
is a unique sliding solution, provided M(r], K) < 1, with mean velocity near threshold given by 



v ~ 



F- F T 



(1 + Vc) 



F - F T 



(30) 



l-M(r),K) ric-r) 

giving (5 = 1, as in the purely elastic case. The critical line rj c (K) separating unique from multi 
valued sliding solutions is determined by M(rj,K) = 1, 

h 2 y 1 

(K + h) 2 / h 



Vc(K) 



1 . 



(31) 



The velocity- force curves and a phase diagram are shown in Fig. |S] for p(h) = e~ h . There is 
a tricritical point at (rj c ,F c = Ft). In contrast to the purely viscous model with K = 0, for 
finite long-time elasticity (K > 0) the behavior is independent of the shape of the pinning force 
distribution, p(h). For r] < rj c , a continuous depinning transition at Ft separates a pinned state 
from a sliding state with unique velocity. In finite dimensions, this transition is likely to remain 
in the same universality class as the depinning of an elastic medium (77 = 0). In our mean-field 
example, the linear response diverges at r/ c , v(rj = r) c ) ~ l/ln(F — Ft). For n > r\ c there is 
hysteresis with coexistence of stuck and sliding states. 
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Numerical simulations of the VE model in two dimensions (dii = d± = 1) indicate a strong 
crossover (possibly a tricritical point) at a critical value of r\ c from continuous to hysteretic depin- 
mnff Although it is always difficult to establish conclusively on the basis of numerics that 

hysteresis survives in the limit of infinite systems, the size of the hysteresis loop evaluated nu- 
merically does appear to saturate to a finite value at large system sizes, indicating that the MF 
approximation may indeed capture the correct finite-dimensional physics. 



V. RELATIONSHIP TO OTHER MODELS AND TO EXPERIMENTS 



Other models of driven systems with inertial-type couplings have been proposed in the literature. 
It is useful to discuss in some detail their relationship to the viscoelastic model considered here. 

In the context of charge density waves, Littlewood ^j] and Levy and collaborators jifil . I~Lt1 | 
modified the Fukuyama-Lee-Rice model Q] that describes the phase of the CDW electrons as an 
overdamped elastic manifold driven through quenched disorder by incorporating the coupling of 
the CDW electrons to normal carriers. This was realized via a global coupling in the equation of 
motion for the phase to the mean velocity of the CDW, not unlike what obtained by a mean-field 
approximation of our viscous coupling. The model was argued to account for the switching and 
non-switching behavior observed in experiments. 

nn 

Schwarz and Fisher |19t 1251 ] recently considered a model of crack propagation in heterogeneous 



solids that incorporates stress overshoot, that is the fact that a moving segment of the crack 
can sometimes overshoot one or more potential static configurations before settling in a new one, 
inducing motion of neighboring segments. These effects may arise from elastic waves that can carry 
stress from one region to another of the driven medium. Stress overshoots, just like topological 
defects in a driven disordered lattice, have an effect similar to that of local inertia and were modeled 
by Fisher and Schwarz by adding couplings to gradients in the local crack velocity in the equation 
of motion for a driven elastic crack. These authors considered an automaton model where time 
is discrete. It is straightforward to define an automaton version of our VE model, where both 
the displacement Uj and time are discrete, as shown in Ref. j^. It is then apparent that the 
automaton version of the viscoeleastic model given in Ref. [3] is identical in its dynamics to the 
model of crack propagation with stress overshoot studied by Schwarz and Fisher, provided the 
strength M of the stress overshoot is identified with the combination 77/(1 + rj). The two models 
differ in the type of pinning considered as the random force used in by Schwarz and Fisher is not 
periodic. We find, however, that the two models have identical mean-field behavior, with a mean- 
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FIG. 6: Mean-field solution of the VE model with a piecewise parabolic pinning potential, p(h) = 8(h — ho) 
and i] = 5. Left frame: phase diagram. Right frame: velocity versus drive for ho/K = 0.5 (blue). ho/K = 1 
(red) and ho/K = 2 (black). Also shown for ho/K = 2 are the discontinuous hysteretic jumps of the velocity 
obtained when F is ramped up and down adiabatically. 

field tricritical point separating continuous from hysteretic dynamical transitions. The connection 
between the viscoelastic and the stress-overshoot model is important because it stresses that distinct 
physical mechanisms (inertia, nonlocal stress propagation, unbound topological defects) at play in 
different physical systems can be described generically by a coarse-grained model that includes a 
coupling to local velocities of the driven manifold. Finally, in a very recent paper, Maimon and 
Schwarz suggested that out of equilibrium a new type of generic hysteresis is possible even when the 
phase transition remains continuous j^H. Driven models with both elastic and dissipative velocity 
couplings may therefore belong to a novel universality class that exhibits features of both first and 
second order equilibrium phase transitions. They clearly deserve further study. 

We now turn briefly to simulations and experiments. For comparison with experiments it is 
useful to point out that the tricritical point of the viscoelastic model can also be obtained by 
tuning the applied force and the disorder strength, rather than the applied force and the viscosity. 
Since the phase diagram does not depend on the form of the disorder distribution, p(h), we choose 
for convenience a sharp distribution, p(h) = 5(h — h®). The phase diagram in the (F,ho) plane 
is shown in Fig. |f)J For weak disorder the depinning is continuous, while for strong disorder it 
becomes hysteretic, with a region of coexistence of pinned and moving degrees of freedom. The 
tricritical point is at (h c ,F c = Ft), with h c = K/ (-^/l + rj — 1). 

Simulations of two-dimensional driven vortex lattices clearly show a crossover as a function 
of disorder strength from an elastic regime to a regime where the dynamics near depinning is 
spatially inhomogeneous and plastic, with coexistence of pinned and moving degrees of freedom 
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39 



40 



52 



531 ] . In fact a bimodal distribution of local velocity was identified in Ref. [40j as 
the signature of plastic depinning. This local plasticity does not, however, lead to hysteresis in 
the macroscopic dc response in two dimension: the mean velocity remains continuous and single- 
valued, although it acquires a characteristic concave-up ward form near depinning that cannot 
be described by the exponent (3 < 1 predicted by elastic models in all dimensions. Hysteresis is, 
however, observed in simulations of three-dimensional layered vortex arrays where the couplings 
across layers are weaker than the in-layer ones jscj • In this case the phase diagram is qualitatively 
similar to that obtained for the viscoelastic model. 

Recent experiments in NbSe2 have argued that memory effects originally attributed in this 



system to "plasticity' of the driven vortex lattice |5J]are actually due edge contamination effects 
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571 ] . In the experiments a metastable disordered vortex phase is injected in a stable 



ordered bulk vortex lattice. Memory effects may then arise in the macroscopic dynamics during 
the annealing of the injected disordered phase. Edge contamination does not, however, explain the 
plasticity seen in simulations, where periodic boundary conditions are used [40j . A possible scenario 
may be that while in the experiments the vortex lattice in the bulk is always in the ordered phase, 
in the simulations the vortex lattice in the bulk of the sample may be strongly disordered even 
in the absence of drive. Such a disordered vortex lattice would then naturally respond plastically 
to an external drive. Finally, it is worth mentioning one experimental situation where hysteresis 
of the type obtained in our model is indeed observed in the macroscopic response. This occurs 
in the context of charge density waves, driven by both a dc and an ac field. In this case the dc 
response exhibits mode- locking steps. The "depinning" from such mode- locked steps was found to 
be hysteretic j^ . 

Several colleagues and students have contributed to various aspects of this work: Alan Middle- 
ton, Bety Rodriguez-Milla, Karl Saunders, and Jen Schwarz. I am also grateful to Jan Meinke for 
help with some of the figures. The work was supported by the National Science Foundation via 
grants DMR-0305407 and DMR-0219292. 



[1] A. J. Liu and S. R. Nagel, Nature 391: 21 (1998); A. J. Liu, this proceedings. 

[2] D. S. Fisher, Phys. Rep. 301: 113 (1998); and references therein. 

[3] G. Blatter et ai, Rev. Mod. Phys. 66: 1125 (1994). 

[4] G. Griiner, Rev. Mod. Phys. 60: 1129 (1988). 



20 



[5] T. Giamarchi and S. Bhattacharya, in High Magnetic Fields: Applications to condensed matter physics 

and spectroscopy, C. Berthier et al., eds. , p. 314 ( Springer- Verlag, 2002). 
[6] E. Y. Andrei, Phys. Rev. Lett. 60: 2765 (1988). 

[7] R. Seshadri and R. M. Westervelt, Phys. Rev. B 46: 5142 & 5150 (1992). 
[8] D. S. Fisher, Phys. Rev. B 31, 1396 (1985). 

[9] T. Giamarchi and P. Lc Doussal, Phys. Rev. Lett. 72: 1530 (1994); Phys. Rev. B 52: 1242 (1995). 
[10] O. Narayan and D. S. Fisher, Phys. Rev. B 46: 11520 (1992). 
[11] F. Nori, Science 271: 1373 (1996). 
[12] A. Tonomura, Micron 30: 479 (1999). 

[13] A. M. Troyanovski, J. Aarts and P. H. Kcs, Nature 399: 665 (1999). 

[14] A. Maeda, et al., Phys. Rev. B 65: 054506 (2002). 

[15] P. B. Littlewood, Solid State Commun. 65: 1347 (1988). 

[16] J. Levy, M. S. Shcrwin, F. F. Abraham, and K. Wiesenfeld, Phys. Rev. Lett. 68: 2968 (1992). 
[17] A. Montakhab, J. M. Carlson, and J. Levy, Phys. Rev. B 50: 11227 (1994). 
[18] M. C. Marchetti, A. A. Middleton, and T. Prellberg, Phys. Rev. Lett. 85: 1104 (2000). 
[19] J. M. Schwarz and D. S. Fisher, Phys. Rev. Lett. 87: 096107 (2001). 

[20] M. C. Marchetti, A. A. Middleton, K. Saunders, and J. M. Schwarz, Phys. Rev. Lett. 91: 107002 
(2003). 

[21] S. H. Strogatz, C. M. Marcus, R. M. Westervelt, and R. E. Mirollo, Phys. Rev. Lett. 61: 2380 (1988); 

Physica D 36: 23 (1989). 
[22] V. M. Vinokur and T. Nattcrmann, Phys. Rev. Lett. 79: 3471 (1997). 

[23] K. Saunders, J. M. Schwarz, M. C. Marchetti, and A. A. Middleton, Phys. Rev. B 69: 37422 (2004). 
[24] M. C. Marchetti and K. A. Dahmen, Phys. Rev. B 66: 214201 (2002). 
[25] J. M. Schwarz and D. S. Fisher, Phys. Rev. E 67: 21603 (2003). 
[26] M. C. Marchetti, Pramana, to be published. 

[27] B. E. Rodriguez-Milla, M. C. Marchetti and A. A. Middleton, unpublished. 
[28] S. Zapperi, J. S. Andrade, Jr. and J. M. Filho, Phys. Rev. B 61, 14791 (2000). 
[29] L. Balents. M. C. Marchetti and L. Radzihovsky, Phys. Rev. B 57: 7705 (1998). 
[30] S. Scheidl, and V. M. Vinokur, Phys. Rev. E 57: 2574 (1998). 
[31] T. Giamarchi and P. Lc Doussal, Phys. Rev. B 57: 11356 (1998). 

[32] Equation ljll|l should include a convective derivative that arises from the transformation to the lab- 
oratory frame [29j. This is negligible near a continuous depinning transition where the mean velocity 
is very small and it drops out in the mean field limit considered here. It will therefore be neglected, 
although it does play a crucial role on the dynamics well into the sliding state [2^] . 

[33] C. R. Myers and J. P. Scthna, Phys. Rev. B 47: 11171 (1993). 

[34] A. A. Middleton, O. Biham, P. B. Littlewood, and P. Sibani: Phys. Rev. Lett. 68, 1586 (1992). 
[35] A. A. Middleton, Phys. Rev. Lett. 68: 670 (1992). 



21 



[36] P. Le Doussal, K. J. Wicse, and P. Chauve, Phys. Rev. E 69: 26112 (2004). 

[37] S. N. Coppersmith, Phys. Rev. Lett. 65: 1044 (1990); Phys. Rev. B 44: 2887 (1991). 

[38] S. N. Coppersmith and A. J. Millis, Phys. Rev. B 44: 7799 (1991). 

[39] A.-C. Shi and A. J. Berlinsky, Phys. Rev. Lett. 67: 1926 (1991). 

[40] M. C. Falcski, M. C. Marchctti and A. A. Middlcton, Phys. Rev. B 54: 12427 (1996). 
[41] L. Balents, M. C. Marchetti, and L. Radzihovsky, Phys. Rev. Lett. 78: 751 (1997). 
[42] A. Zippclius, B. I. Halperin, and D. R. Nelson, Phys. Rev. B 22: 2514 (1980). 

[43] See, for instance, D. R. Nelson in Phase Transitions and Critical Phenomena, C. Domb and J. Lcbowitz 

eds. (Academic, London, 1983), Vol. 7, pp. 76-79. 
[44] M. C. Marchctti and K. Saunders, Phys. Rev. B 66: 224113 (2002). 
[45] J. P. Boon and S. Yip, Molecular Hydrodynamics (McGraw-Hill, New York, 1980). 
[46] The strain-stress relation is nonlocal in both space and time. Only the nonlocality in time is retained 

here. 

[47] A. E. Koshelev and V. M. Vinokur, Phys. Rev. Lett. 73, 3580 (1994). 

[48] K. Dahmcn and J.P. Sethna, Phys. Rev. Lett. 71: 3222 (1993); Phys. Rev. B 53: 14872 (1996). 

[49] R. da Silveira and M. Kardar, Phys. Rev. E 59: 1355 (1999). 

[50] C. J. Olson, C. Reichhardt and V. M. Vinokur, Phys. Rev. B 64: 140502 (2001). 

[51] R. Maimon and J. M. Schwarz, Phys. Rev. Lett. 92: 255502 (2004). 

[52] C. J. Olson, C. Reichhardt, and F. Nori, Phys. Rev. Lett. 81: 3757 (1998). 

[53] N. Gronbcch-Jensen, A. R. Bishop, and D. Dominguez, Phys. Rev. Lett. 76: 2985 (1996). 

[54] S. Bhattacharya and M. J. Higgins, Phys. Rev. Lett. 70: 2617 (1993); M. J. Higgins and S. Bhat- 

tacharya, Physica C 257: 232 (1996). 
[55] Y. Paltiel, E. Zeldov, Y. N. Myasoedov, K. Shtrikman, S. Bhattacharya, M. J. Higgins, Z. L. Xiao, E. 

Andrei, P. L. Gammcl, and D. J. Bishop, Nature 403: 398 (2000). 
[56] Y. Paltiel, Y. Myasoedov, E. Zeldov, G. Jung, M.L. Rappaport, D. E. Feldman, M.J. Higgins, and S. 

Bhattacharya, Phys. Rev. B 66: 60503 (2002). 
[57] M. Marchcvsky, M. J. Higgins, and S. Bhattacharya, Phys. Rev. Lett. 88: 87002 (2002). 
[58] M. J. Higgins, A. A. Middleton, and S. Bhattacharya Phys. Rev. Lett. 70: 3784-3787 (1993) 



